Extremal property of a simple cycle 
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We study systems with finite number of states Ai (i = 1, . . . , n), which obey the first order kinetics 
(master equation) without detailed balance. For any nonzero complex eigenvalue A we prove the 
inequality j^j < cot — . This bound is sharp and it becomes an equality for an eigenvalue of a 
simple irreversible cycle Ai — > A2 — > . . . — > A n — > Ai with equal rate constants of all transitions. 
Therefore, the simple cycle with the equal rate constants has the slowest decay of the oscillations 
among all first order kinetic systems with the same number of states. 

PACS numbers: 82.20.Rp, 82.40.Bj, 02.50.Ga 



Damped oscillations in kinetic systems. We study 
systems with finite number of states Ai (i = 1, . . . , n). A 
non- negative variable Pi corresponds to every state Ai. 
This pi may be probability, or population, or concentra- 
tion of Ai. For the vector with coordinates (pi) we use 
the notation P. 

We assume that the time evolution of P obey first order 
kinetics (the master equation) 



(i) 



where qij (i,j = 1, ...,n, i ^ j) are non-negative. In 
this notation, is the rate constant for the transition 
Aj — > Ai. Any set of non- negative coefficients (i ^ j) 
corresponds to a master equation. 

Let us rewrite ((T|) in the vector notations, P = KP, 
where the matrix K = (fey) has the elements 



ki 



Qij 



if i = j. 



(2) 



The standard simplex A„ = {P \ Pi > 0, ^2iPi = 1} 
is forward-invariant with respect to the master equation 
([T]) because it preserves positivity and has the conserva- 
tion law "YliiVi = const. This means that any solution 
of (p} P(t) with the initial conditions P(to) G A„ re- 
mains in A n after to- Pit) € A n for t > to- One can use 
this forward invariance to prove some known important 
properties of K . For example, there exists a non-negative 
vector P* G A„ [p* > 0) such that KP* =0 (equilib- 
rium). Indeed, any continuous map of A„ — > A„ has a 
fixed point, therefore exp(Kt) has a fixed point in A„ 
for any t > 0. If exp(Kt)P* = P* for some P* E A„ 
and sufficiently small t > 0, then KP* — because 
exp{Kt)P = P + tKP + o{t). 

The proof that K has no nonzero imaginary eigenval- 
ues for the systems with positive equilibria gives another 
simple example. We exclude one zero eigenvector and 
consider K on the invariant hyperplane where J2iPi ~ 0- 



If K has a nonzero imaginary eigenvalue A, then there 
exists a 2D if-invariant subspace L, where K has two 
conjugated imaginary eigenvalues, A and A = —A. Re- 
striction of exp(ivt) on L is one-parametric group of ro- 
tations. For the positive equilibrium P* the intersection 
(L + P*) n A„ is a convex polygon. It is forward invari- 
ant with respect to the master equation (JTJ because L is 
invariant, P* is equilibrium and A„ is forward invariant. 
A polygon on a plane cannot be invariant with respect to 
one-parametric semigroup of rotations exjp(Kt) (t > 0). 
This contradiction proves the absence of imaginary eigen- 
values. We use the reasoning based on forward invariance 
of A„ below in the proof of the main result. 

The master equation obeys the principle of detailed 
balance if there exists a positive equilibrium P* (p* > 0) 
such that for each pair i,j (i ^ j) 



HjPi = QjiPi ■ 



(3) 



After Onsager [T|, it is well known that for the systems 
with detailed balance the eigenvalues of K are real be- 
cause K is under conditions ([3]) a self-adjoined matrix 
with respect to the entropic inner product 



(see, for example, @, Q). 

Detailed balance is a well known consequence of mi- 
croreversibility. In 1872 it was introduced by Boltzmann 
for collisions [4(. In 1901 Wegscheider proposed this prin- 
ciple for chemical kinetics [5]. Einstein had used this 
principle for the quantum theory of light emission and 
absorbtion (1916, 1917). The backgrounds of detailed 
balance had been analyzed by Tolman Q. This princi- 
le was studied further and generalized by many authors 



Systems without detailed balance appear in applica- 
tions rather often. Usually, they represent a subsystem 
of a larger systems, where concentrations of some of the 
components are considered as constant. For example, the 
simple cycle 
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Ai A 2 



A.. 



Ai 



(4) 



2 



is a typical subsystem of a catalytic reaction (a catalytic 
cycle). The complete reaction may have the form 



S + At -> A 2 



P, 



(•5) 



on this plane has the following matrix 



JC 



5RA 



3A 9?A 



where 5 is a substrate and P is a product of reaction. 

The irreversible cycle (j3| cannot appear as a limit of 
systems with detailed balance when some of the constants 
tend to zero, whereas the whole catalytic reaction §E§ 
can The simple cycle (01 can be produced from the 
whole reaction ([SJ if we assume that concentrations of S 
and P are constant. This is possible in an open system, 
where we continually add the substrate and remove the 
product. Another situation when such an approximation 
makes sense is a significant excess of substrate in the sys- 
tem, [S] ^> [Ai] (here we use the square brackets for the 
amount of the component in the system). Such an excess 
implies separation of time and the system of intermedi- 
ates {Ai} relaxes much faster than the concentration of 
substrate changes. 

In systems without detailed balance the damped oscil- 
lations are possible. For example, let all the reaction rate 
constants in the simple cycle be equal, qj+ij = qin = q > 
0. Then the characteristic equation for K is 

det(K - XI) = (-q - A) n + q n (-l) n+1 = 

and A = — q+qexp (^77^) (k — 0, . . . ,n— 1). For nonzero 
A, the ratio of imaginary and real parts of A is 



In this coordinate system 



2-rrfc I 



I5RAI 



1 — COS 



2irk I 



nk 

cot 

n 



7T 

< cot — . 

n 



The maximal value, cot — , corresponds to k = 1. For 
large n, cot — m — and oscillations in the simple cycle 
decay rather slowly. 

Estimate of eigenvalues. Let us consider the general 
master equation ([T]) without any assumption of detailed 
balance. 

Theorem 1. For every nonzero eigenvalue A of matrix 
K 



I3AI 7T 

W\- COt n 



(6) 



Proof. Let us assume that the master equation ([TJ has a 
positive equilibrium P* G A„: for alii = 1, . . . , n p* > 
and 



£ top} = Yl 



QjiPi 



The systems with non-negative equilibria may be consid- 
ered as limits of the systems with positive equilibria. 

Let A be a complex eigenvalue of K and let L be a 
2D real subspace of the hyperplane J2i Pi — that cor- 
responds to the pair of complex conjugated eigenvalues, 
(A, A). Let us select a coordinate system in the plane 
L + P* with the origin at P* such that restriction of K 



exp(f/C) 



exp(iSA) cos(^A) - exp(iKA) sin(*3A) 
exp(m\) sin(^A) exp(t3?A) cos(iSA) 



The intersection A = (L + P*) n A„ is a polygon. It 
has not more than n sides because A„ has n (n — 2)- 
dimensional faces (each of them is given in A„ by an 
equation pi =0). For the transversal intersections (the 
generic case) this is obvious. Non-generic situations can 
be obtained as limits of generic cases when the subspace 
L tends to a non-generic position. This limit of a se- 
quence of polygons cannot have more than n sides if the 
number of sides for every polygon in the sequence does 
nor exceed n. 

Let the polygon A have m vertices Vj (m < n). 
We move the origin to P* and enumerate these vectors 
Xj = Vj — P* anticlockwise (Fig [TJ . Each pair of vectors 
Xj,Xj_|_i (and x m ,xi) form a triangle with the angles on, 
Pi and 7i, where is the angle between Xj and Xj+i, and 
f3 m is the angle between x m and xi. The Sine theorem 

|xi I |xi+i I Ix™ I Ixi I 

gives J " — 1 + 1 "'' — 1 11 

sin n 



sin ■ji 



sin 71 



Several elementary identities and inequalities hold: 

< Cii,/3i,ji < 7r; y j f3 i = 2ir; ai + f3 l + j t = ir; 

i 

jPJsin<x; = J^Jsin7i (the closeness condition). 



(7) 



These conditions ([7]) are necessary and sufficient for the 
existence of a polygon A with these angles which is star- 
shaped with respect to the origin. 

Let us consider the anticlockwise rotation (Q'A < 0, 
Fig. H]). The case of clockwise rotations differs only in 
notations. For the angle S between K~x.i and Xj, sini5 = 
— 5A, cos 8 = — 5RA and tan S = i£. 

For each point x 6 L + P* (x ^ P*), the straight 
line {x + e/Cx |e G 1} divides the plane L + P* in two 
half-plane (Fig.[lJ dotted line). Direct calculation shows 
that the semi-trajectory {exp(£/C)x 1 1 > 0} belongs to 
the same half-plane as the origin P* does. Therefore, if 
6 < ai for alH = 1, . . . , m then the polygon A is forward- 
invariant with respect to the semigroup exp(i/C) (t > 0). 
If 8 > ai for some i then for sufficiently small t > 
exp(£/C)xi ^ A because /Cx^ is the tangent vector to the 
semi-trajectory at t = 0. 

Thus, for the anticlockwise rotation (3A < 0), the 
polygon A is forward-invariant with respect to the semi- 
group exp(</C) {t > 0) if and only if 8 < ai for all 
i = 1, . . . , m. The maximal 8 for which A is still forward- 
invariant is <5 max = minijaj}. We have to find the poly- 
gon with m < n and the maximal value of miiii{ai}. Let 
us prove that this is a regular polygon with n sides. Let 



3 




FIG. 1: The polygon A is presented as a sequence of vectors 
x;. The angle /3j between vectors x, and x;+i and the angles 
on and ji of the triangle with sides Xj and Xj+i are shown. 
In the Fig., rotation goes anticlockwise, i.e. SJA < 0. In this 
case, the polygon A is invariant with respect to the semigroup 
exp(t/C) (t > 0) if and only if 8 < on for all i = 1, . . . , m, 
where 5 is the angle between the vector field /Cx and the 
radius- vector x. 



us find the maximizers ai,Pi,~fi (i = l,...,m) for the 
optimization problem: 

minjcvi} — » max subject to conditions ([7]). (8) 

i 

For solution of this problem, all at are equal. To prove 
this equality, let us mention that mhii{ai} < ? under 
conditions Q (if all o» > 5 then the polygonal chain A 
cannot be closed). Let min^ ai = a. Let us substitute 
in ([7]) the variables which take this minimal value by 
a. The derivative of the left hand part of the last condi- 
tion in ([Jj with respect to a is not zero because a < j. 
Assume that there are some ctj > a. Let us fix the val- 
ues of fli (i — 1, . . . , m). Then 7, is a function of ai, 
7i = 7r — Pi — ai. We can use the implicit function the- 
orem to increase a by a sufficiently small number e > 
and to change the non-minimal ay by a small number 
too, aj i — y ctj — 0', = 0(e). Therefore, at the solution of 
([5]) all aj — a (j — 1,. . . , m). 

Now, let us prove that for solution of the problem (|SJ) 
all ^ are equal. We exclude 7$ from conditions ([TJ and 
write Pi + a < ir; < (3i, a; 



Let us consider this equality as equation with respect to 
unknown a. The function log sin a; is strictly concave on 
(0,7r). Therefore, for X{ € (0,n) 



log sin 



1 



E 



Xi > 



1 



E 



log sin Xi 



and the equality here is possible only if all equal. 
Let a* 6 (0, 7r/2) be a solution of (|9]). If not all the values 
of Pi are equal and we replace Pi in © by the average 
value, P = 2£ then the value of the right hand part of (|9]) 
increases and sin a* < sin(/3 + a*). If we take all the Pi 
equal then ([9]) transforms into elementary trigonometric 
equation sin a = sin(/3 + a). The solution a of equation 
(J9]) increases when we replace Pi by the average value: 
a > a* because sin a* < sin(/3 + a*), a G (0, 7r/2) and 
sin a monotonically increases on this interval. So, for the 
maximizers of the conditional optimization problem ((5J) 
all /3j = — and ai = 7i = 5 — 5-. The maximum of a 
corresponds to the maximum of m. Therefore, m = n. 
Finally, max{J} = | — - and 



f |9fA| 1 7T 

max fe) =cot ^ 

This is exactly the same value as for an eigenvalue of the 
simple cycle of the lengths n with equal rate constants, 
A=-«(l-exp(^)). □ 



m log sin a 



logsin(/3 4 + a). 



(9) 



Discussion. The simple cycle with the equal rate con- 
stants gives the slowest decay of oscillations or, in some 
sense, the slowest relaxation among all first order kinetic 
systems with the same number of components. The ex- 
tremal properties of the simple cycle with equal constants 
were noticed in numerical experiments 25 years ago [9[. 
V.I. Bykov formulated the hypothesis that this system 
has extremal spectral properties. This paper gives the 
answer: yes, it has. 
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